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The rise of the total, elastic and inelastic hadronic cross sections at high energies is investigated 
by means of an analytical parametrization, with the exponent of the leading logarithm contribution 
as a free fit parameter. Using derivative dispersion relations with one subtraction, two different fits 
to proton-proton and antiproton-proton total cross section and rho parameter data are developed, 
reproducing well the experimental information in the energy region 5 GeV - 7 TeV. The parametriza- 
tion for the total cross sections is then extended to fit the elastic (integrated) cross section data in 
the same energy region, with satisfactory results. From these empirical results we extract the energy 
dependence of several physical quantities: inelastic cross section, ratios elastic/total, inelastic /total 
cross sections, ratio total-cross-section/elastic-slope, elastic slope and optical point. All data, fitted 
and predicted, are quite well described. We conclude that, once the TOTEM Collaboration results 
for the total and elastic cross sections at 7 TeV are confirmed, our analysis indicates: (1) an in- 
crease of the hadronic cross sections with the energy faster than the log-squared bound by Froissart 
and Martin; (2) asymptotic limits 1/3 and 2/3 for the ratios elastic/total and inelastic/total cross 
sections, respectively, a result in agreement with unitarity. These indications corroborate recent 
theoretical arguments by Ya. I. Azimov on the rise of the total cross section. 



PACS: 13.85.-t Hadron-induced high- and super-high-energy interactions, ll.10.Jj Asymptotic prob- 
lems and properties, 13.85.Lg Total cross sections. 
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I. INTRODUCTION 

The rise of the hadron-hadron total, elastic and inelastic cross sections at high energies is an ex- 
perimental fact. However, the exact energy dependence involved and a widely accepted theoretical 
explanation for this increase have been open problems for a long time. Presently, these aspects con- 
stitute one of the greatest challenges for QCD, the gauge field theory of the strong interactions. The 
total cross section is connected with the forward elastic amplitude (optical theorem) and therefore, 
perturbative techniques cannot access the large distances phenomena that permeate the elastic chan- 
nel. On the other hand, non-perturbative approaches are not yet able to evaluate soft scattering states 
from first principles. As a consequence, the dependence of the total hadronic cross section with the 
energy can not yet be directly predicted by QCD. 

In the phenomenological context, a wide variety of models present good descriptions of the available 
data. However, since they are characterized by distinct physical pictures [lHf|, not exclusively related 
to QCD, a global, unified and well accepted approach is still missing. 

In the theoretical context, formal results from axiomatic local quantum field theory, general princi- 
ples and further arguments, state that the total cross section cannot grow with the energy faster than 
log-squared (ThU, a result also extended to the inelastic cross section (Iol - [r^ |. However, as recently 
discussed by Azimov [l3l - fl5j , it is not obvious if these derivations can be directly applied to hadronic 
processes (QCD) and formal arguments suggest that the total cross section could grow faster than 
log-squared. Moreover, recent result from lattice QCD indicates an universal asymptotic logarithm- 
squared dependence for the hadronic total cross section, once some specific assumptions are made 
[lq . That is, this lattice QCD result does not represent an unique or exclusive solution. 

In the experimental context, new results from the LHC and new estimates from the Pierre Auger 
Observatory for the proton-proton (pp) total cross sections are expected to shed light on the subject. 
In particular, the first result for the pp total cross section at 7 TeV (LHC) has been obtained by the 
TOTEM Collaboration, through a luminosity- dependent measurement [171 ]. Although, in principle, 
an experimental result considered rather conservative, due to agreement with some model predictions, 
further investigations brought out some unexpected aspects. In fact, once included in amplitude anal- 
ysis this point can not be described by standard parametrization with leading log-squared dependence 
on the energy, as first discussed by us in [l8| (hereafter referred to as FMS) and recently (2012) also 
verified by the Particle Data Group fl9| . 

The main ingredient in the FMS analysis has been the use of an analytical parametrization intro- 
duced by Amaldi et al. in the seventies [2(| and also used by the UA4/2 Collaboration, in the nineties 
[2lj . This parametrization is characterized by power contributions at low energies (reggeons) and a 
leading contribution at high energies (pomeron), parametrized by a power law in the logarithm of the 
energy with the exponent as a free real parameter. Detailed analysis, with different methods and fit 
variants to total cross section data only, from pp and antiproton-proton (pp) scattering above 5 GeV, 
led to the conclusion that the real exponent in the logarithmic term is statistically consistent with 2 if 
the TOTEM datum is not included in the analysis, but above 2 (around 2.2 - 2.3) if this point takes 
part in the fitted data-set. 

In this work we extend the FMS analysis on pp and pp scattering above 5 GeV in several aspects, 
all of them including the result of the TOTEM Collaboration at 7 TeV: 

1. As discussed in [18], since the analytical parametrization is nonlinear in several parameters, the 
data reductions demand start (feedback) values for the fit parameters [22j and depending on 
the choices the fit result may be different. Here, two different choices have been considered, one 
from the FMS analysis [HI and another one based on the recent results from the PDG [T3| ■ 

2. Using derivative dispersion relations, we connect the total cross section, o~tot, with the p param- 
eter (ratio between the real and imaginary parts of the forward amplitude). With this analytical 
formalism, we first present the predictions for p(s) from individual fits to a tot data and after 
that we develop novel simultaneous fits to Otot and p data. In both cases (individual and simul- 
taneous fits) we discuss in some detail the important role of the subtraction constant, embodied 
in the dispersion relations. 

3. As commented above, the log-squared bound, originally stated for the total cross section by 
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Froissart-Martin-Lukaszuk [7h1| , 



00 < 4 1- 2 f- 



) 



where s is the center-of-mass energy squared, s r a constant, the pion mass, has been recently 
extended to the inelastic cross section by Martin- Wu- Roy-Singh |10l - [l3 |. 




Therefore, from s-channel unitarity, a similar bound is expected for the elastic (integrated) cross 
section, er e ;(s). Based on the connection between the total cross section and the forward elastic 
scattering amplitude (optical theorem) and since the FMS analysis of the total cross section has 
indicated a power larger than 2, we address here the extension of the same parametrization to 
the elastic cross section. Namely, with start (feedback) values of the simultaneous fits to o~tot 
and p data we develop new empirical fits to a e i data, with the same analytical parametrization. 

4. With the results from the empirical fits for cr tot (s), p(s) and o~ e i(s), we present predictions for 
the energy dependence of several physical quantities: ai ne i, the ratios o- e i/o~tot, &inel/o~tou the 
optical point, as well as a result for the ratio between o~tot and the elastic slope B and from that 
for the proper slope B. 

5. Here, all physical quantities that are compared with experimental data (fitted and predicted 
quantities), are displayed with uncertainty regions, evaluated by means of standard (analytical) 
error propagation from the free fit parameters [22| . 

All the experimental data considered (fitted and predicted) are quite well described. We conclude 
that, if further measurements of the pp total cross section confirm the present result from the TOTEM 
Collaboration at 7 TeV, our analysis leads to solutions that favor an increase of the total and elastic 
cross sections faster than the log-squared bound. In that case, the violation of the bound does not 
imply in violation of unitarity, as recently discussed by Ya. I. Azimov. Indeed, we obtain that, 
asymptotically, cr e i/o- tot — > 1/3 and Omei/o'tot — > 2/3. However, although statistically consistent and 
corresponding to quite good descriptions of the existing data, our results do not constitute unique 
solutions, but possible solutions for the rise of the hadronic cross sections. 

The manuscript is organized as follows. In Sect. II we treat the fits to o~tot and p data: the analytical 
parametrization for a tot is introduced and the analytical connection with p by means of dispersion 
relations is presented and discussed, with emphasis on the role of the subtraction constant. In Sect. 
Ill the extension of the parametrization to elastic cross sections data is presented. In Sect. IV we 
display the predictions for the inelastic cross section, the ratios involving cross sections, the optical 
point, the elastic slope and compare with the corresponding experimental data, followed by discussions 
on each obtained result. Our conclusions and some critical remarks are the contents of Sect. V. 



At high energies, in terms of the elastic scattering amplitude F, the total cross section and the p 
parameter in the forward direction can be expressed by 





II. TOTAL CROSS SECTION AND THE RHO PARAMETER 



0- t ot{s) 



lm.F(s,t = 0) 



(Optical Theorem) 



(1) 



s 



Pis) 



ReF(s,t = 0) 
Im F(s,t = 0)' 



(2) 



where t is the four momentum transfer squared. 
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A. Experimental Data 

Our interest here is the investigation of the high and asymptotic energy region, associated with 
particle-particle and antiparticle-particle scattering. For that reason we shall consider only elastic 
collisions with the highest energy interval in terms of available data, namely pp and pp scattering. 
With this restrictive choice we do not take into account data from other reactions, available only 
in the regions of intermediated and low energies or any constraint dictated by a supposed universal 
behavior. 

The input data set for fits concerns only accelerator data, covering the region from 5 GeV up 
to 7 TeV. However, estimations of the pp total cross section from cosmic-ray experiments will be 
displayed as illustrative results. The reason why we do not include this information in the fits is the 
model-dependence involved in the extraction of the proton-proton total cross section from proton-air 
production cross sections 123 -261. The accelerator data on a t ot and p from pp and pp scattering have 
been compiled fro m Il7|. 1271 l28j and the estimations of the pp total cross sections from cosmic-ray 
experiments from j28l . |29| . 



B. Analytical Parametrization 

The analytical parametrization to be used here in the fits to pp and pp total cross sections data 
(and in Sect. Ill to elastic cross sections data), is given by (l8l [20l |2~D| 



a pp (s) 
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(4) 



where si = 1 GeV 2 (fixed), ax, b\, a 2 , b 2 (low energies), a, /?, 7 and Sh (high energies) are, in general, 
free real fit parameters. 

For the case of total cross section all the parameters above have specific physical interpretations 
in the context of the Regge-Gribov theory. They are associated with reggeon exchanges (first two 
terms) and pomeron exchanges (last two terms) at low and high energy regions, respectively (l8l. [30|. 
In the case of the elastic cross section we consider this parametrization as an empirical ansatz or as 
a representation. 



C. Dispersion Relations and the Subtraction Constant 



Amplitude analyses of the total cross section rising usually include the information on the p pa- 
rameter through fits to both quantities [3fJ-|39|. The connection between atot and p is obtained by 
means of dispersion relations (integral and/or derivatives forms), or the asymptotic prescriptions for 
crossing even and odd amplitudes (Phragmen-Lindeloff theorem) [40L l4lj | . In the high-energy region, 
of interest in this work, the dispersion relations demand one subtraction [42l |43| , which means the 
addition of one more fr ee p arameter in the simultaneous investigation of both quantities. 

In the FMS analysis [Tsj] the fits have been developed only with the o~tot data, without the inclusion 
of the p information. Here, using derivative dispersion relations with one subtraction, we also treat 
simultaneous fits to both quantities. To this end, in this subsection, we shortly review some aspects 
related to integral and derivative dispersion relations, with emphasis on the origin and role of the 
subtraction constant. In the next subsections, we present our predictions for p{s) from fits to atot 
data and after that, we develop and discuss new simultaneous fits to p and a to t data. 
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Integral and Derivative Dispersion Relations 



Following the notation in [26| , for pp and pp scattering, analyticity and crossing symmetry allow us 
to connect <Jtot{s) and p(s) through two compact and symmetric formulas: 



p pp [s)a 



E(a + ) + 0(a-), 



(5) 



^{s)of ot {s) = E{cx + )-0{a-l 
where the even (+) and odd (— ) cross section are related to the physical cross sections by 



(6) 



a±(s) 



i'i> 



(7) 



and E(a+), 0{a ) arc analytic transforms connecting the real and imaginary parts of crossing even 
and odd scattering amplitudes, respectively (associated with cr + and These analyticity relations 
are usually expressed in an integral form (Hilbert transform (44|) and in the case of the forward 
direction the standard once-subtracted integral dispersion relations (IDR), with poles removed, may 
be expressed by [4ll - l43l l45j | 



K 2s f°° 
E int (a+) = - + — P ds' 



(8) 



p oo 


' s' ' 




s /2 _ s 2 





(9) 



where K is the subtraction constant and P denotes principal value. We recall that in the even case, 
IDR with one and two subtractions are equal and in the odd case, those without subtraction and 
with one subtraction are also equal 4l|. Therefore, with our parametrization of interest (3-4), the 



subtraction in the even part is adequate since, asymptotically, from equation (7), we have: 



In 7 (s) and 



as 



We shall return to the subject of the subtraction constant in the next subsection. 

On the other hand-for classes of functions of interest, integral dispersion relations can be replaced 
by derivative forms }46l448| | . known as derivative dispersion relations (DDR), which may be more 
useful for some practical calculations. In these formulas, differentiation with respect to the logarithm 
of the energy occurs in the argument of a trigonometric operator, as in the standard form deduce d by 
Bronzan, Kane, and Sukhatme in the high-energy approximation (so — > in Eqs. (8) and (9)) |49l.l5C 



E d 
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Oder(c-) = tan 
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(11) 



corresponding, as before, to singly-subtracted DDR. 

This completes the analytical approach: with an input parametrization for o~f% t (s), a PP t (s), Eqs. 
(5), (6) and (7) allow, in principle, the determination of pP p (s) and p pp (s) by means of either the IDR, 
(8) and (9), or the DDR, (10) and (11) in the high-energy approximation. 

In our case, with the In 7 (s/sh) term and 7 a free real parameter, IDR demand numerical methods 
and therefore a nonanalytic approach for error propagation from the fit parameters. The use of 
prescriptions seems to us unjustified in the region of intermediate and low energy data (5 - 20 GcV), 
since they are asymptotic results 40]. On the other hand, DDR allow an analytical approach and even 
the high-energy results, (10) and (11), can be applied in the low energy region, an effect related to 
the presence of the subtraction constant K . Since this effect plays a fundamental role in our analysis, 
let us first recall the origin, meaning and role of this constant, with focus on its disadvantages and 
advantages in the present case. 
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2. The Role of the Subtraction Constant 



Subtractions are useful tools to deal with the convergence of IDR and they are essentially based on 
a change of variable. In general grounds, suppose a given dispersion relation for a complex function 
f{z) and let Xq be some point on the real axis, at which f{z) is analytic. The corresponding relation 
for [f(z) — f(xo)]/[z — xo] leads to the subtracted dispersion relation 0, H3, [5l| . Separating the real 
an imaginary parts of f(z), the IDR for Kef(z) includes now the constant Re/(a;o), the subtraction 
constant. Therefore, in a particular physical problem it is necessary to determine Re/(xo) for some 
xo- For scattering amplitudes xq — > s = is assumed, so that for an even amplitude or cross 

section as in Eq. (8), the subtraction constant reads 

K = ReF+(s = 0), 



meaning the analytic continuation of the even amplitude at s = 4l| . 

Therefore, the subtraction constant has a well justified mathematical origin, related to the con- 
vergence of the integral. On the other hand, we cannot attribute to it a (real) physical meaning or 
value, since it is associated with the nonphysical region. However, in applying dispersion relations to 
connect a to t and p this constant appears explicitly and therefore must be taken into account. Note 
that to assume K — is only a particular (and unjustified) choice, as any other. In this respect two 
disadvantages are in order to be recalled: 

1. As an unknown constant, it plays the role of one more free parameter. However, the lack of a 
physical meaning contrasts with all the other free parameters present in the analytical parametriza- 
tion of the total cross section, namely reggeons intercepts, reggeons strengths and leading pomeron 
contribution. 

2. In data reductions, as a free parameter, its value affects the final values of the others physical 
free parameters. Moreover, it may seem that, once appearing in the form K/s in Eqs. (8) and (10), 
its influence is limited to the low energy region. However, that is not true: typical (global) fits are 
characterized by strong correlation among all free parameters. That means K may affect the fit results 
even at high and asymptotic energies, as already discussed and demonstrated in [26l IHol. l52j| . 

Summarizing, once introduced in the fit procedure, one may, in principle, select and even anchor 
the outputs. That was one of the reasons why simultaneous fits, including the p information, were 
not considered in the FMS analysis (see however Appendix A in j!8j). 

On the other hand, in applying derivative dispersion relations deduced in the high-energy approx- 
imation (previous subsection), the subtraction constant develops an useful property, related to the 
equivalence between DDR and IDR. These aspects have been already demonstrated and discussed in 
pa . [HoL [52j and in what follows we recall the main points. 

As commented before, the DDR deduced by Bronzan-Kane-Sukhatme (10-11) have been obtained 
in the high-energy approximation, which means considering sq — > in the lower limit of the IDR 
(8) and (9). Since for pp and pp scattering, s = 4m^ ~ 3.5 GeV 2 , it could be expected that, for 
a given parametrization, the results obtained by means of IDR with so = 4m 2 fixed and through 
DDR are not the same. That, however, is not always the case. In [5,0] and references therein, the 
practical equivalence of IDR and DDR has been investigated in the case of the Donnachie-Landshoff 
parametrization, with both degenerate [53| and non-degenerate [H, H3| meson trajectories. Several 
variants have been considered in simultaneous fits to a to t and p data from pp and pp scattering: (i) 
energy cutoffs at 5 and 10 GeV; (ii) subtraction constant K = and if as a free fit parameter; (iii) 
DDR and IDR with both so = Am^ and so — 0. In all cases, once the subtraction constant is used as 
a free parameter, the numerical results obtained through DDR (high-energy approximation) and IDR 
with fixed sq = 4m? (without the high-energy approximation) are the same, up to 3 significant figures 
in the fit parameters and x 2 /DOF. This effect is a consequence of the absorption of the high-energy 
approximation in an effective subtraction constant. Specifically, as analytically demonstrated in Sect. 
4.4 of [j5(|, the coefficient A of the series expansion associated with the integration from s' = to 
s' = 4mp is absorbed by the subtraction constant, since the series can be put in the form 

[K + A] - + O (\ 
s \s z 

This effect of absorption of the high-energy approximation by the subtraction constant has also been 
confirmed by other authors [55|,[56|]. Therefore, despite the disadvantages of the subtraction constant, 
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it can have here a useful practical or pragmatic role as a free fit parameter in the DDR: it leads to the 
same result that could be obtained through the IDR without the high-energy approximation. 

We recall that, for classes of functions of interest in high-energy elastic scattering, DDR can be 
analytically extended down to 4 - 5 GeV [H, H3] or even below and for the whole energy interval 
(above the physical threshold), either in the form of a double infinity series, as first deduced by Avila 
and Menon [Ha . |59| . or in the form of a single series, as demonstrated by Ferreira and Sesma [60| . 
However it seems to us simpler here to consider the pragmatic role of the subtraction constant, since 
as explained above, with this constant as a free parameter the results obtained with the DDR (10) 
and (11) (high-energy approximation) are the same as those obtained through the IDR (8) and (9) 
(without the high-energy approximation). Therefore, let us focus on the DDR (10) and (11) with K 
as a free fit parameter. 



3. Analytical results with Derivative Dispersion Relations 



The applicability of derivative dispersion relations in amplitude analysis has been critically reviewed 
by Avila and Menon in special the replacement of the trigonometric operators in (10) and (11) 
by the corresponding series (see also [Hj]). In practice, with the presence of the ln 7 (s/s/j) term, the 
derivative transforms (10) and (11) can be evaluated through the operational expansion, introduced 
by Kang and Nicolescu (6l| (see also Sect. 3 in (50l|): 



Eder(o~+) — 



K 



1 



2 dins 3 V 2 dins 



d 



15 V 2 dins 



0+0), 



(12) 



Oder (O'- 




er- (s) > dhis. 



(13) 



With parametrization (3-4) as input and from (7), the evaluation of the power terms (low energy 
contribution) results in closed forms (sum of the series) . For the logarithm term (leading contribution) 
a fast convergence is already obtained up to third order, since as we shall show 7 ~ 2.2 — 2.4. In this 
case we obtain 18] 



, . K 
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where 



A = -M, 
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/ 3 7 [ 7 _1][ 7 _ 2 ], 



C = l5 



/3 7 [ 7 -l][ 7 -2][ 7 -3][ 7 -4] 



(15) 



and for the odd contribution, 



-62 

(16) 

With the above results, (14,15,16) and parametrization (3-4) for af^ t and af^ t , Eqs. (5) and (6) 
lead to the analytical expressions for p pp (s) and pP p (s): 



Oder(cr-) = -0,2 COt 
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with ^4, B and C given by Eq. (15). Note that the analytical results imply that ass4 oo: 

1 



p oc 



In s 



0, 



(17) 



(18) 



which is in agreement with the rigorous asymptotic result by Khuri and Kinoshita [62j |. 



D. Fits and Results 



In this subsection, after discussing general aspects of the fit procedures and methodology, we present 
the results, through Eqs. (17) and (18), of individual fits to atot data with corresponding predictions 
for the p{s) and those from simultaneous fits to Otot and p(s) data. A detailed discussion of all the 
results here obtained is presented in the next subsection. 



1. Fit Procedures and Methodology 

The nonlinearity of the fit (constants b\, 62, 7 and Sh in Eqs. (3-4)) demands a choice of the initial 
(feedback) values of the free parameters [22]. Here we shall address two different choices, denoted and 
explained in what follows. 

Fit 1. In the FMS analysis, based only on the total cross section data, two ensembles, two methods 
and six variants have been tested for different feedback values and different fixed and free 
parameters. Here we shall select as a representative result with 7 > 2 the one denoted in [l8[ as 
Ensemble Vs max = 7 TeV, Method 2 and Variant 5 (Table IV and Figure 6 in pj). The reason 
for this choice is five-fold: 

(i) x 2 / D F closest to 1; 

(ii) smaller relative error in the exponent 7; 

(iii) result closest to the TOTEM point at 7 TeV; 

(iv) the numerical value of the exponent 7 does not correspond to the largest or the smallest 
value obtained in the analysis (Figure 8 in [3). 

(v) the reggeon intercepts, 61 and 62, are fixed to 0.5, so that we can investigate the effects of 
this assumption (compared with Fit 2 discussed below, in which the intercepts are not fixed). 

Fit 2. The 2012 PDG analysis [l^] is based on the hi ghes t rank parametrization result by the COM- 
PETE Collaboration, introduced ten years ago |30l . l3~fl | . The functional form can be seen as a 
particular case of parametrization (3-4) in which the exponent 7 is fixed to 2. The analysis by 
PDG includes different reactions and tests on the universatility. Since this work corresponds to 
another fit result, with a kind of "conservative character" (namely 7 = 2), we have considered 
the values of the parameters they have obtained as start values for fits with parametrization 
(3-4). We note that in doing so we are applying the parametrization to a subset of the data 
analyzed by the PDG and by the COMPETE Collaboration. 
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The data reductions have been performed with the objects of the class TMinuit of ROOT Framework 
[63j j and all results correspond only to the cases of full convergence. As tests of goodness of fit we 
consider the x 2 per degrees of freedom (DOF), with the x 2 directly obtained from the MINUIT Code. 
In all fits we have adopted the Confidence Level (CL) of rj 68 % (one standard deviation), which 
means that the projection of the \ 2 distribution in (TV + l)-dimensional space {N — number of free 
fit parameters) contains 68 % of probability (22j . 

In the MINUIT Code, the correlation matrix gives a measure of the correlation between each pair of 
free parameters, with numerical limits ± 1 (full correlation) and (no correlation) [22] (we shall return 
to this point in what follows). The error matrix provides the variances and covariances associated 
with each free parameter. This information is used in the analytic evaluation of the uncertainty 
regions in the fitted and extracted (predicted) physical quantities through standard error propagation 
procedures [22T |. 

In what follows we shall address some of the points raised on the role of the subtraction constant in 
subsection II. C. 2, now with the results related to parametrization (3-4), Fits 1 and 2 and the derivative 
dispersion results (17) and (18). First we treat the predictions for p(s) from fits to er tot data and after 
that, simultaneous fits to both <r tot and p data. As explained, the consistent applicability of the 
DDR demands the subtraction constant if as a free fit parameter. However, in order to illustrate 
some aspects of this constant we shall consider two variants in all cases: if = and if as a free fit 
parameter. 

2. Individual Fits to Total Cross Section Data and Predictions for the Rho Parameter 

• Total Cross Section 

In the case of Fit 1, the values of the parameters obtained in 18], through parametrization (3-4) 
and fit to only the total cross section data, are displayed in the second column of Table I (atot data 
column) , with the corresponding statistical information on the data reduction. The dependencies of the 
cross sections on energy are shown in Figure 1, including uncertainty regions from error propagation, 
together with the experimental data analyzed (accelerator) and estimations of a PP t from cosmic-ray 
experiments. 

In the case of Fit 2 to total cross section data, the feedback values of the parameters and the fit 
results are displayed in the 5th and 6th columns of Table I, respectively (a to t data columns). In this 
case, the parameter sm that appears associated with the power contributions in the PDG analysis, 
has been included in the reggeon strengths a\ and a-i- We distinguish however the scale factor in the 
logarithm contribution (sh)- The results for the total cross section, including uncertainty regions, and 
experimental information are shown in Figure 2. 

• Predictions for p(s) 

From Eqs. (17) and (18), any prediction for p pp (s) and p pp (s) from the total cross section fit 
depends on the value of the subtraction constant. As commented before, we shall consider the two 
variants if = and if as a free fit parameter to the p data. Specifically, we use as input the fixed 
values of the parameters for the total cross section (Table I) , with si — 1 GeV 2 . For if — (fixed) we 
obtain the direct prediction for p(s). In the case of if as a free parameter we use if = as start value 
and fit only the p data with this parameter, by fixing all the other parameters (from the total cross 
sections). In this case the first run with the MINUIT Code allows us to obtain the corresponding 
X 2 /DOF for if = 0. This statistical information and the values of if as fit parameter are displayed in 
the 3rd and 4th columns of Table I with Fit 1 and 7th and 8th column of Table I with Fit 2 (p data 
columns). The dependencies of p on energy, including uncertainty regions from error propagation and 
the experimental data, are shown in Figure 1 with Fit 1 and Figure 2 with Fit 2 (both for if = and 
if as a free fit parameter). 

3. Simultaneous Fits to Total Cross Section and Rho Data 

For simultaneous fits to a to t and p data, we have used as feedback the values of the parameters 
from Fit 1 and Fit 2 to a tot-, displayed in Table I (2nd and 6th columns). Here, we also consider the 
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two variants: K = and K as a free fit parameter (in this case, with K = as start value). The fit 
results and statistical information are displayed in Table II. The results for <7tot(s) and p(s) from Fit 
1 are shown in Figures 3 (K — 0) and 4 [K as free parameter) and from Fit 2 in Figures 5 (K = 0) 
and 6 (K as free parameter). 

TABLE I: Results of the individual fits to u to t data and the predictions for the p parameter, with the subtrac- 
tion constant K = and K as a free parameters and Fits 1 and 2 (see text). Also included are the statistical 
information on the data reductions: degrees of freedom (DOF) and reduced x 2 . The parameters ai, 0,2, a 
and /3 are in mb, sn in GeV 2 and bi, 62, 7 and K are dimensionless. In all cases, si = 1.0 GeV 2 (fixed) m 
parametrization (3-4). 





Fit 1: 






Fit 2: 








Ensemble: 


a to t data 


p data 




atot data 




p data 






[18] 


K =0 


K free 


Initial values 
[19] 


fit results 


K =0 


K free 


Oi 


56.5±1.1 






46.06 


60.2±1.5 






61 


0.5 (fixed) 






0.462 


0.4850±0.0079 






a 2 


27.70±0.28 






34.0 


33.1 ± 1.7 






62 


0.5 (fixed) 






0.550 


0.540 ± 0.015 






a 


33.65±0.22 






34.71 


31.52 ± 0.22 






P 


0.1301±0.0086 






0.265 


0.0575± 0.0025 






7 


2.213±0.024 






2.0 


2.422 ± 0.016 






Sh 


3.90±0.52 






16.2 


0.566 ± 0.065 






K 







44.8±4.9 









38.5±4.9 


DOF 


158 


76 


75 




156 


76 


75 


x 2 /DOF 


0.967 


2.67 


1.58 




0.934 


2.28 


1.48 



TABLE II: Results of the simultaneous fits to a to t and p data, with feedback values from Fit 1 and Fit 2 
(Table I), for K = and if as a free fit parameter. Same legend as Table I. 





From Fit 1: 




From Fit 2: 






K = 


K free 


K = 


K free 


ai 


48.58 ± 0.72 


51.5 ± 4.3 


65.99 ± 0.55 


59.8 ± 1.3 


61 


0.5 (fixed) 


0.5 (fixed) 


0.2758 ± 0.0026 


0.4541 ± 0.0069 


Ct2 


26.97 ± 0.27 


27.63 ± 0.28 


34.5 ± 1.5 


34.1 ± 1.6 


b 2 


0.5 (fixed) 


0.5 (fixed) 


0.5511 ± 0.0094 


0.547 ± 0.010 


a 


35.41 ± 0.15 


35.16 ± 0.89 


5.59 ± 0.24 


29.78 ± 0.22 


P 


0.264 ± 0.041 


0.32 ± 0.11 


0.2575 ± 0.0028 


0.0693 ± 0.0025 


7 


2.018 ± 0.057 


1.952 ± 0.098 


1.8769 ± 0.0040 


2.346 ± 0.013 


Sh 


18.7 ± 3.8 


22 ± 15 


0.00903 ± 0.00058 


0.383 ± 0.041 


K 





41. 2± 6.5 





33.5 ± 5.9 


DOF 


232 


231 


232 


231 


X 2 /DOF 


1.37 


1.15 


1.14 


1.11 



E. Discussion and Partial Conclusions 



Beyond the distinct start values for the parameters in Fits 1 and 2, an essential difference concerns 
the reggeon intercepts, which are fixed to 61 = 62 = 0.5 in the former case and free in the latter. 
Based on this feature and the results presented in Tables I and II and Figures 1 to 6, we have the 
comments that follow, including some partial conclusions. 
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• Individual Fits to Total Cross Section Data and Predictions for the Rho Parameter (Table 1, Figures 
1 and 2) 

1. Both fits indicate 7 > 2 with good statistical confidence (Table I): 
Fit 1: 7 - 2.21 ± 0.02 and x 2 /DOF - 0.97, 

Fit 2: 7 - 2.42 ± 0.02 and x 2 /DOF - 0.93. 

2. The uncertainty regions for the total cross sections from both fits are in agreement with the 
TOTEM result at 7 TeV (Figures 1 and 2). 

3. In the predictions for the p parameter, the subtraction constant affects the results in the low 
energy region, y^s < 20 GeV, essentially a change of curvature in pP p around 10 GeV, in the 
case of K free (Figures 1 and 2). 

4. In all cases, the global description of the p data is quite good and, as expected, the best statistical 
results are obtained with K as a free fit parameter (Table I). 

5. All the results with Fit 2 corroborate those obtained in the FMS analysis. 

• Simultaneous Fits to Total Cross Section and Rho Data (Table II and Figures 3 to 6) 

1. Compared with individual fits to a to t data, the inclusion of p data in simultaneous fits leads to 
a decrease in the value of the exponent 7 in all cases analyzed (Tables I and II): 

Fit 1: 7 ~ 2.21 => 7 ~ 2.01 (K = 0) and 7 - 1.95 (K free) 

Fit 2: 7 - 2.42 7 ~ 1.88 (K = 0) and 7 - 2.35 (K free) 

2. The results with Fit 1 are consistent with the log-squared bound (Table II): 
7 - 2.02 ± 0.06 for K = and 7 - 1.95 ± 0.10 for K free. 

The uncertainty regions for the total cross sections barely include the TOTEM result at 7 TeV. 
(Figures 3 and 4). 

3. The result with Fit 2 for K = is consistent with the log-squared bound, 7 <~ 1.877 ± 0.0004, 
but the uncertainty region does not describe the TOTEM data at 7 TeV (Figure 5). 

4. The result with Fit 2 for K free is consistent with an increase of the total cross section faster 
than the log-squared bound, 7 ~ 2.35 ± 0.01 and the uncertainty region describes the TOTEM 
data at 7 TeV (Figure 6). 

As commented before, we have included the results with K = only to illustrate some aspects of the 
subtraction constant, outlined above. Since the formal and consistent applicability of the DDR here 
used demands the subtraction constant as a free fit parameter, we shall focus our partial conclusions 
on this case only. 

Concerning Fit 1, we understand that fixing the reggeon intercepts to 0.5, results in a substantial 
constraint in the possible increase of the total cross section, when the p data is included in simultaneous 
fits. That is a consequence of the correlations among all the fit parameters, in special between K and 
7. On the other hand, with Fit 2, the reggeon intercepts and strengths are free to control the low 
and intermediate energy data, leaving the parameters associated with the constant and the leading 
logarithmic contributions adequately applied to the highest energy data. The correlation matrix for 
both fits with if as a free fit parameter is displayed in Appendix A, where the correlation between 
K and 7 is stressed. Moreover, the best statistical result has been obtained with Fit 2 (x 2 /DOF 
closest to 1) and within the uncertainty regions all the experimental information on <j tot and p is well 
described (Figure 6). 

Based on these facts, we shall select as our representative result that one obtained with Fit 2 in 
the simultaneous fits to a to t and p data in the case of if as a free fit parameter (5th column in Table 
II). This result, as well as those obtained with individual fits to a to t data (Fits 1 and 2) support an 
increase of the total cross section faster than the log-squared bound. 
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III. FITS TO ELASTIC CROSS SECTION DATA 

From s-channel unitarity, an increase of the total cross section faster than log-squared is expected 
to be, in some way, connected with a similar effect in the elastic and/or inelastic cross sections. At 
least that is what is suggested by the log-squared bound of the total cross section, recently extended 
to the inelastic cross section, as recalled in our introduction. Therefore, to go one step further, a point 
of interest is to test the applicability of the same parametrization (3-4) to the elastic or inelastic cross 
section data. In order to address this possibility we shall consider only the elastic cross section for 
the four reasons that follow. 

(i) From the optical theorem the total cross section is connected with the elastic amplitude in the 
forward direction. 

(ii) Total and elastic cross sections can be directly measured or estimated without the need of model 
assumptions. That however is not the case for the inelastic cross section due to the single and 
double dissociative contributions. 

(iii) The least ambiguous way to estimate the inelastic cross section is from unitarity, namely o~i ne i = 
Ctot ~ o~ei [64j j . which suggests a more fundamental character for the last two quantities. 

(iv) To select a statistically consistent ensemble of <Ji ne i data for fit is not an easy task due to both 
model-dependence in "direct" estimations and the different data from different experiments and 
energies for <r tot and o~ e i , if using unitarity. 

The experimental data on the elastic cross section from pp and pp scattering have been compiled 
from [28j. For the fits to the elastic cross section data with parametrization (3-4) we have used as 
start values those obtained with our representative result: Fit 2 and K as a free fit parameter in the 
simultaneous fits to o~tot and p data. 

In this case, a crucial point concerns the exponent 7 that controls the asymptotic behavior of the 
cross sections. Given that we have obtained for the total cross section 7 tot = 2.346 ± 0.013, a value 
different from this for the elastic cross section, j e i, has no physical meaning, due to direct violation of 
unitarity by the ratio o- e i/o~tot or the reduction of this ratio to zero at asymptotic energies. Therefore 
we shall explore the possibility that the same exponent applies for the elastic cross section as well. 
Specifically, we fix this parameter to the central value obtained for the total cross section, j e i = 2.346. 
With the above assumptions and fit procedures, we obtain the results displayed in Table III and 
Figure 7, corresponding therefore to a good reproduction of the experimental data. We shall discuss 
this result together with predictions for other physical quantities in the next Section. 

TABLE III: Results of the fit to a e i data through parametrization (3-4), with feedback values from Fit 2, K 
free (5th column in Table II). Same legend as Table I. 



ai 


30.7 ± 3.6 


bi 


0.551 ± 0.037 


ai 


0.236 ± 0.071 


b 2 


0.134 ± 0.012 


a 


4.28 ± 0.14 


P 


0.02358 ± 0.00054 


7 


2.346 (fixed) 


Sh 


0.978 ± 0.023 


DOF 


97 


X 2 /DOF 


1.62 



IV. PREDICTIONS FOR OTHER PHYSICAL QUANTITIES 

Based on the empirical fits developed for atot(s), p(s) and a e i(s), we present here the predictions 
for several physical quantities, followed by a discussion on each result. 
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A. Inelastic Cross Section and Ratios Elastic/Total and Inelastic/Total 

From unitarity, by subtracting parametrization (3-4) with the corresponding values of the param- 
eters displayed in Table II for atot(s) and Table III for <r e /(s), we obtain our predictions for ai„ e i(s). 
The result is shown in Figure 8 together with the experimental data [28], [65l - l67l | and the uncertainty 
region. With this result and those for atot(s) and a e i(s), we extract the corresponding predictions for 
the ratios a e i/atot and <Ji ne i/otot, displayed in Figures 9 and 10, respectively. 



B. Optical Point, Elastic Slope and the Ratio with the Total Cross Section 

In this section we first recall the definitions of the elastic slope B, the optical point and an expression 
connecting the ratios <J e i/&tot and <T tot /B. After that we present our predictions. 
In terms of the elastic amplitude, the differential cross section is expressed by 



da 



(s,t) 



1 



\F(s,t)f 



dt""'"' 167TS 2 

and the slope of the elastic differential cross section in the forward direction is defined as 



" d 






(**)] 


di 





B(s,t = 0) 



From (1), (2) and (19), the optical point is given by 

dt t _ n 16-7T 



t=0 



The integrated elastic cross section reads 

CT e z(s) 



da 

, dt 



(s, t)dt, 



(19) 



(20) 



(21) 



(22) 



where to defines the physical (kinematic) region. 

Concerning the differential cross section, experimental data indicate a sharp forward peak, fol- 
lowed by a dip-bump or dip-shoulder structure above ~ 0.5 GeV 2 (Tevatron, LHC). Typically, these 
structures are located more than 5 decades below the optical point, Eq. (21). These experimental 
facts are important in the determination of the integrated elastic cross section, since in this case 
the differential cross section can effectively be represented by an exponential fall off, simulated by a 
model- independent parametrization [17} . 



da 
~dt 



da 
~dt 



-Bl 



t=0 



with B the (constant) forward slope. In that case, by assuming to 
clastic cross section is given by the approximate result 



(23) 

-co in Eq. (22), the integrated 



a el (s) = -4t ^M[l + p 2 ] 
V ' B(s) 16tt 1 H 1 

and therefore, 

o-tot(s) 167r a e i(s) 



B{s) [1 + p 2 ] at t(s)' 
which is very close to the MacDowell-Martin bound [g8| , 



(24) 



(Ttotjs) a e i(s) , , 

— -— < 18?r (25) 

B(s) <Jtot{s) 
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a result recently discussed in more detail in |69j. For our purposes, the main ingredient in the result 
(24) is the possibility to investigate the behavior of atot{s)/ B(s) from the information on the ratio 
Oei{s)/otot{s) (see also [69|]). 

With our empirical results for er tot (s), p{s) and u e i(s) we can predict the optical point (21), the 
ratio atot/B (24) and from that the elastic slope B. The results for the last two quantities are shown 
in Figures 11 and 12, respectively. The predictions for several quantities of interest at 7, 8, 14 and 57 
TeV are displayed in Table IV, together with the TOTEM results at 7 TeV. 



TABLE IV: TOTEM results at 7 TeV and predictions for several physical quantities associated with the pp 
elastic scattering at the LHC and Pierre Auger Observatory energies. 



Physical 


7 TeV 


7 TeV 


8 TeV 


14 TeV 


57 TeV 


quantity 


TOTEM [HI 










Otot (mb) 


98.3 ± 2.8 


96.40 ± 0.97 


98.7 ± 1.0 


108.6 ± 1.2 


137.0 ± 1.9 


P 




0.1362 ± 0.0016 


0.1357 ± 0.0016 


0.1333 ± 0.0015 


0.1261 ± 0.0013 


a e i (mb) 


24.8± 1.2 


24.31 ± 0.31 


25.03 ± 0.33 


28.18 ± 0.39 


37.26 ± 0.59 


Uinei (mb) 


73.5 ± 3.1 


72.1 ± 1.0 


73.7 ± 1.1 


80.4 ± 1.3 


99.7 ± 2.0 


O e t/o t ot 


0.252 ± 0.014 


0.2522 ± 0.0041 


0.2536 ± 0.0042 


0.2595 ± 0.0046 


0.2720 ± 0.0057 


Oinel/otot 


0.7477 ± 0.038 


0.7478 ± 0.0041 


0.7465 ± 0.0042 


0.7405 ± 0.0046 


0.7280 ± 0.0057 


Otot IB 


12.56 ± 0.43 


12.44 ± 0.21 


12.51 ± 0.21 


12.82 ± 0.23 


13.46 ± 0.28 


B (GeV" 2 ) 


20.10 ± 0.36 


19.89 ± 0.39 


20.25 ± 0.40 


21.76 ± 0.46 


26.15 ± 0.66 


do/dt\ t =o (mbGeV" 2 ) 


504 ± 27 


483.6 ± 9.7 


507 ± 10 


613 ± 14 


974 ± 27 



C. Discussion 



From Figure 8 the prediction for Oi ne i(s) presents good agreement with all the accelerator data and 
barely describes the recent Auger estimation at 57 TeV. At 7 TeV it favors the TOTEM data, a result 
obtained through unitarity and therefore corresponding to a model-independent evaluation. 

The results for the ratios of elastic/total and inelastic/total cross sections are displayed in Figures 
9 and 10. Within the uncertainty region, the predictions are also consistent with the TOTEM results 
at 7 TeV. These results lead to the consequences that follows. 

Asymptotically (s — » oo), from parametrization (3-4) and by denoting the parameters associated 
with atot and a e i fits by the corresponding indexes, we have 



Qei _^ Pel 
Otot Plot 



and 



Oinel 
Otot 



From Tables II and III we obtain 



Oex_ 

Otot 



-» 0.340 ±0.015 



which are consistent with the rational limits: 

o e i 1 
► q 

^tot o 



and 



and 



Oinel 
Otot 



Oinel 
Otot 



Pel 
Ptot 



0.660 ±0.015, 



(26) 



indicating, therefore, a constant asymptotic value for the ratio o e i/atot which lies below 1/2, the 
black-disk limit. 

It is interesting to note that the above results corroborate recent arguments attributing the black- 
disk limit to a combination of the soft processes, namely elastic and diffractive. Specifically, denoting 
Odiff the single and double diffractive contributions, Grau-Pacetti-Pancheri-Srivastava have proposed 
the following asymptotic limit [701 ] 



Pel 
Otot 



a diff 
Otot 



as 



oo. 
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This argument is based on the formulation of eikonal models [6J, l71| and according to our results it 
is expected: 

a diff _^ 1 
<Jtot 6 ' 

Another aspect concerns the recently proposed empirical parametrization for the ratio elastic/total 
cross section by Fagundes and Menon [691 ] . 

— (s) = Atanh( 7 i + 72 lns + 7 3 ln 2 s), (27) 

where 7i, i = 1, 2, 3 are free fit parameters and A represents the asymptotic limit. In this work 
the authors have considered two extrema cases, A = 1/2 (black-disk limit) and A = 1 (beyond the 
black-disk limit), with good descriptions of the pp data above 10 GeV. Our rational limit (26) implies 
A = 1/3 and therefore another scenario, lying asymptotically below the black-disk limit, as already 
pointed out. 

The predictions for the ratio utot/B through Eq. (24) and for the elastic slope B are displayed in 
Figures 11 and 12, respectively. We recall that, even in the forward direction (t = 0), the determination 
of the slope from the differential cross section data demands some interval in the momentum transfer. 
That can explain some of the discrepant points in both figures since the corresponding intervals from 
different experiments are not always the same. Within the uncertainty regions, the predictions show 
good agreement with the experimental information at high energies, in special with the TOTEM result 
at 7 TeV. Asymptotically, from Eq. (24) and since p — > (Sect. II. C. 3), our rational limit (26) gives 

Otot 167T 

and therefore, an asymptotic increase of the slope also faster than log-squared. As commented before, 
this result represents a different scenario from that discussed in [6^| and also in [72j . 

In Table IV we display our numerical results (with uncertainties) for all the physical quantities here 
discussed and at different energies of interest in present and future experiments. At 7 TeV, except for 
the p parameter not determined by the TOTEM Collaboration, the compatibility of our results with 
the experimental data is quite good, for all physical quantities. 

At last, another aspect deserving some comments concerns the predictions for the proton-proton 
total cross sections beyond the (accelerator) energy region analyzed, namely the estimations of a^ t 
from proton-air production cross section in cosmic-ray experiments. As already commented, these 
points did not take part in our data reductions and have been displayed only as qualitative illustrations. 
The estimations at the highest energies (Fly's Eye and Auger Collaborations) depend on extrapolation 
from phenomenological models tested only in the accelerator energy region. This fact, associated with 
the relative small cosmic-ray flux above ~ 20 TeV, results in extremely large uncertainties, as shown in 
our figures. Even though, the fit extrapolation suggests good agreement with these points (Figure 6), in 
special our numerical prediction at 57 TeV (Table IV) with the recent Auger Collaboration estimation 
for the total cross section, a^ t = [133 ± 13 (stat)[^ 20 (sys) ± 16 (Glauber)] mb [29j . However, it is 
important to note that systematic and theoretical (Glauber) uncertainties are equally likely quantities, 
that is, do not follow Gaussian distributions, as is the case of statistical uncertainties. That means 
the above central value is equally likely to lye in any place limited by the corresponding non-statistical 
uncertainties. For example, added in quadrature the resulting uncertainty corresponds to ~ 23 mb 
(around 17 % of the central value). This result contrasts with the small systematic uncertainty of 
the TOTEM point a 7 TeV, namely 2.8 mb (around 2.8 %). Therefore, despite the suggestive global 
agreement, the large equally likely uncertainty may imply in different scenarios, as already discussed 
in [17]. We add that this observation puts limits on recent arguments by Block and Halzen in what 
concerns the description of the Auger Collaboration result quoted above [77J (see also our critical 
remarks at the end of the next section) . 



V. CONCLUSIONS AND CRITICAL REMARKS 



In this work an amplitude analysis on the forward pp and pp elastic scattering in the energy region 5 
GeV - 7 TeV has been presented. The main point concerns the investigation of the log-squared bound 



D. A. Fagundes, M. J. Menon, P. V. R. G. Silva 



16 



in the rise of the cross sections with the energy. To address this question, an analytical parametrization 
introduced by Amaldi et al. (2pj . with the exponent of the leading logarithm contribution as a 
free fit parameter, has been used. Analytical connection between a tot and p parameter has been 
obtained by means of subtracted derivative dispersion relations. The presentation has been divided 
in three parts: (1) individual and simultaneous fits to a tot and p data, with focus on the role of the 
subtraction constant; (2) extension of the parametrization to the elastic cross section data with fixed 
exponent of the logarithm contribution; (3) predictions for other physical quantities and study of some 
consequences. 

As we have shown and discussed in detail all analyzed and predicted data are quite well described. 
The results indicate: (i) an increase of the hadronic cross sections faster than the log-squared bound; 
(ii) asymptotic limits for the ratios elastic/total and inelastic/total cross sections consistent with 1/3 
and 2/3, respectively. Therefore, the violation of the bound does not imply in violation of unitarity, 
as recently discussed by Ya. I. Azimov: the fast rise may correspond to a rapid high-energy increase 
of the scattering amplitude in non- physical regions [H, [3 ■ 

The possible violation of the log-squared bound stated by Froissart and Martin, may look like as 
an unexpected breakdown of a paradigm. However, that is not the case in what concerns the use of 
parametrization (3-4), namely the possibility to treat 7 as a free parameter and not to fix it to 2 (or 
1). In fact, up to our knowledge, before the FMS analysis, this parametrization was used in two works 
only, the first one by Amaldi et al. in the seventies and the second by the UA4/2 Collaboration in 
the nineties. Both analyses treated only pp and pp elastic scattering and simultaneous fits to atot and 
p data. In the energy interval 5 GeV< ^/s < 62.5 GeV, Amaldi et al. obtained [2(| 

7 = 2.10 ± 0.10 

and the UA4 Collaboration, with the data set extended up to 546 GeV [2lj . 

~ - 25+ - 35 

Therefore, both results, well known for a long time, indicate a rise of the total cross section faster 
than the log-squared bound. 

In the FMS analysis two ensembles on a to t data have been tested, one with data up to \fs max — 
1.8 TeV and another one including the TOTEM result, namely \fs max — 7 TeV. The novel aspect of 
the analysis consisted in the evidence that with the former ensemble the 7 value is consistent with 2, 
but above 2 with the latter one. In the case of y/s max — 7 TeV, several solutions have been obtained 
with different methods and variants 181: 



7 sa 2.10 ± 0.03 (method 1, variant 1), 7 « 2.27 ± 0.04 (method 1, variant 3), 

7 rj 2.21 ± 0.02 (method 2, variant 5), 7 « 2.10 ± 0.11 (method 2, variant 6). 

Although the fits had been restricted to atot data only, the use of derivative dispersion relations led 
to predictions for the p parameter consistent with the experimental information, even with the largest 
7 values [la ]. Here, we have obtained with Fit 2 to atot data 

7 « 2.42 ±0.02 

and in the case of simultaneous fits to atot and p data, 

7 w 2.35 ±0.01. 

Therefore, based on all the above results, we understand as convincing evidence the fact that 
parametrization (3-4) applied to experimental data may lead to statistically consistent solutions with 
7 > 2, possibly in the interval 2.2 - 2.4. 

The leading contribution in parametrizations of the total cross section is expected to be associated 
with pomeron exchanges [HQ- For 7 = 1, the constant plus Ins terms correspond to a double pole 
at J = 1 and for 7 = 2 a triple pole (expressed by In 2 s, Ins and the constant terms) [30|. 45 1. The 



case of non-integer exponent, with < 7 < 2 corresponds to a strong-coupling scenario 73|, |74| and 



a fractional power, 7 = 3/2 (in general 1 < 7 < 2), is phenomenologically indicated by the eikonal 
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minijet model with soft gluon £; t -resummation [64], [7(1 [75], [76( . Presently, we do not have a physical 
interpretation for a leading logarithm contribution with the exponent 7 as a real parameter greater 
than 2. Perhaps thinking about an effective exponent, representing the sum of different contributions 
might not be so speculative. However, to discuss this point may yet be premature for the critical 
comments that follows. 

Once associated with the leading contribution in parametrization (3-4), the 7 value from data 
reductions is essentially determined by the experimental information at the highest energy interval 
available. With y/s = 1.8 TeV we are faced with the well known discrepant values of dtot obtained 
by the CDF, E710 and E811 Collaborations. This set of data gives no practical information on 
the rise of the total cross section, except for a possible effective mean value in data reductions. It 
seems peculiar the fact that most phenomenological models and even amplitude analyses present 
consistence with this imaginary mean value. On the other hand, compared with all the experimental 
data available, the TOTEM results at 7 TeV are characterized by extremely small uncertainties, as 
can be easily seen in the figures. Our analysis has been essentially based on fits to total cross section 
and elastic cross section, for which the uncertainties in the TOTEM data are around 2.8 % and 4.8 
%, respectively. Therefore, in our data reductions, with the reduced x 2 as a test of goodness of fit, 
these two points have played an essential role in the determination of the 7 parameter. In other 
words, we understand that our evidence for 7 > 2 is a consequence of the TOTEM results at 7 TeV, 
as also demonstrated in 18]. The 2012 PDG analysis corroborates this conclusion since, as already 
commented in our introduction, the COMPETE log-squared parametrization does not describe the 
TOTEM data (Figure 46.10 in 

Therefore, based on the above critical comments, we understand that confirmation of the TOTEM 
results at 7 TeV is necessary before any firm conclusion can be drawn on the 7 value and consequently 
on the results here presented and discussed. In this sense, the expected measurement of the total cross 
section at 8 TeV through a luminosity-independent method seems to be crucial. We also emphasize 
that our results represent possible consistent statistical solutions for the behavior of the total cross 
section, but do not correspond to unique solutions. Despite these strong limitations, we hope this 
analysis and our previous work [18[ might contribute with further interpretations and developments in 
the investigation of the rise of the hadronic cross sections at high energies, a problem, in our opinion, 
still unsolved. 
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Appendix A: Correlation Matrices 



TABLE V: Correlation Matrix [HH from Fit 1: 7 free parameters, 61 = b 2 = 0.5 (fixed). 



Parameter 


0,1 


ai 


a 


P 


7 


Sh 


K 


Ol 


1.000 


0.063 


-0.990 


-0.786 


0.696 


-0.929 


0.558 


a 2 


0.063 


1.000 


-0.009 


0.063 


-0.073 


0.027 


0.271 


a 


-0.990 


-0.009 


1.000 


0.852 


-0.770 


0.968 


-0.510 


P 


-0.786 


0.063 


0.852 


1.000 


-0.987 


0.954 


-0.304 


7 


0.696 


-0.073 


-0.770 


-0.987 


1.000 


-0.899 


0.247 


Sh 


-0.929 


0.027 


0.968 


0.954 


-0.899 


1.000 


-0.424 


K 


0.558 


0.271 


-0.510 


-0.304 


0.247 


-0.424 


1.000 
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TABLE VI: Correlation Matrix [H, HI from Fit 2: 9 free parameters. 



Parameter 


Ol 


bi 


<Z2 


&2 


Q 


P 


7 


Sh 


K 


0,1 


1.000 


0.787 


0.208 


0.176 


-0.077 


0.061 


-0.139 


-0.252 


0.382 


bi 


0.787 


1.000 


0.026 


0.020 


0.457 


-0.062 


-0.131 


-0.129 


0.317 


0.2 


0.208 


0.026 


1.000 


0.976 


-0.063 


-0.005 


0.041 


0.063 


-0.206 


b 2 


0.176 


0.020 


0.976 


1.000 


-0.059 


-0.008 


0.039 


0.062 


-0.266 


a 


-0.077 


0.457 


-0.063 


-0.059 


1.000 


-0.139 


0.178 


0.557 


0.029 




0.061 


-0.062 


-0.005 


-0.008 


-0.139 


1.000 


-0.816 


0.364 


0.014 


7 


-0.139 


-0.131 


0.041 


0.039 


0.178 


-0.816 


1.000 


0.120 


-0.031 


Sh 


-0.252 


-0.129 


0.063 


0.062 


0.557 


0.364 


0.120 


1.000 


-0.057 


K 


0.382 


0.317 


-0.206 


-0.266 


0.029 


0.014 


-0.031 


-0.057 


1.000 
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FIG. 1: Individual fit to a to t data through parametrization (3-4) from Fit 1 and predictions for p(s) with K 
= and K as a free fit parameter (Table I). 
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FIG. 2: Individual fit to a to t data through parametrization (3-4) from Fit 2 and predictions for p(s) with K 
= and K as a free fit parameter (Table I). 
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FIG. 3: Simultaneous fits to a to t and p data from Fit 1 and K — (Table II). 
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FIG. 5: Simultaneous fits to a to t and p data from Fit 2 and K — (Table II). 
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FIG. 11: Predictions for the ratio between total cross section and the slope parameter and experimental data. 
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FIG. 12: Predictions for the slope parameter and experimental data. 



